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qj , Abstract. We study analytically, in one dimension, the survival probability P s (t) 

up to time t of an immobile target surrounded by mutually noninteracting traps 
I ■ each performing a continuous-time random walk (CTRW) in continuous space. We 

consider a general CTRW with symmetric and continuous (but otherwise arbitrary) 
c/3 ' jump length distribution f(r]) and arbitrary waiting time distribution ip( T )- The traps 

+-* . are initially distributed uniformly in space with density p. We prove an exact relation, 

£H valid for all time t, between P s (t) and the expected maximum E\M(t)] of the trap 

'T 1 ■ process up to time t, for rather general stochastic motion x tlap (t) of each trap. When 

*^ , £trap(£) represents a general CTRW with arbitrary f{rf) and i/j(t), we are able to 

compute exactly the first two leading terms in the asymptotic behavior of E[M(t)] 
f\ i for large t. This allows us subsequently to compute the precise asymptotic behavior, 

P s (t) ~ a exp[— bt e ], for large t, with exact expressions for the stretching exponent 9 

fvj ' and the constants a and b for arbitrary CTRW. By choosing appropriate f(rf) and - 0( r ): 

^ . we recover the previously known results for diffusive and subdiffusive traps. However, 

O^ | our result is more general and includes, in particular, the superdiffusive traps as well 

j-^s ■ as totally anomalous traps. 
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1. Introduction 

The study of the survival probability of a target particle surrounded by diffusing traps 
is a classic problem in physical chemistry that goes back to Smoluchowski in 1916 
PQ. Smoluchowski introduced this problem to compute the rate of chemical reactions 
in the diffusion-driven regime. The basic model and its variants [21 E] have since 
found a large number of applications in a wide variety of contexts such as in reaction- 
diffusion systems [3J, chemical kinetics [HI El [3 El El [101 EH], predator-prey models in 
population dynamics [121 E] and also in a wide class of the so called 'walker persistence' 
problems [T3] . 

In the simplest setting of the model known as the 'target annihilation' problem [TJ] , 
the target particle is immobile at the origin and the noninteracting traps, initially 
distributed uniformly over space with a finite density p, each perform an unbiased 
continuous-time Brownian motion with a diffusion constant D. When any of these 
traps hits the target particle at the origin, the target particle gets annihilated. The 
average survival probability P s (t) of the target up to time t, where the average is over 
the initial positions of the traps distributed uniformly over space with density p, has been 
computed exactly [H] in all dimensions. In particular, in the simplest one dimensional 
case, P s (t) decays with t as a stretched exponential for all t, 



P s (t) = exp 




(1) 



In d = 2, P s (t) ~ exp[— pc2t/lnt] while for d > 2, P s (t) ~ exp[— pcdt] for large t, where 
the constants C2 and Cd are known [TJ]. In d > 2, the target particle needs to have a finite 
size while for d < 2, the target can be considered to be a point particle. In the original 
Smoluchowski version of the problem, the quantity F(t) = — ln[P s (i)] has a physical 
interpretation: it represents the net flux of traps to a static absorber at the origin up 
to time t 0EI [ID]. In d = 1, the net flux F(t) = ^pJ^ thus grows algebraically ~ t 1/2 
with time t. 

The original Smoluchowski or equivalently the 'target annihilation' problem has 
been generalized in the literature in several ways, e.g., (a) when the target itself moves 
stochastically (not necessarily diffusive) while the noninteracting traps perform normal 
diffusion and (b) when the target is static but the traps, while still noninteracting, 
perform, in general, non-diffusive stochastic motion (c) when both the target and the 
traps move nondiffusively in general and (d) the famous Donsker-Varadhan problem [15] 
where the target moves diffusively but the traps are immobile. In this paper, we will 
focus on the cases (a) and (b) only. 

An interesting variant belonging to class (a) above is the 'diffusive target 
annihilation' problem where the target particle itself diffuses with a diffusion constant 
D [16J. This problem has seen a recent flurry of activity [HI HH1 [191 EDI EI] and it 
has been proved rigorously that in one dimension, the average survival probability P s (t) 
of the diffusing target particle has the same asymptotic stretched exponential decay 
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as in Eq. (QQ) (and hence is independent of D ), though with considerable subleading 
corrections to the leading behavior for intermediate t. Another variant in class (a) is the 
'ballistic target annihilation' problem where the target particle moves ballistically with a 
constant velocity c. In this case, the survival probability decays faster, P s (t) ~ exp[— M] 
in all dimensions d, where the inverse decay rate $ has been computed exactly for d < 2 
andd = 3 [22]. 

The variants in class (b) where the target is static and the traps undergo subdiffusive 
motion or (c) where both the target and the traps undergo subdiffusive motion have 
been studied extensively by Yuste and collaborators [231 1211 [25j I26J [27] . 

Another recent extension belonging to class (b) above is to the case where the 
target is static, each trap performs independent Brownian motion but resets to its own 
initial position with a constant rate r [28J. This problem has interesting implication in 
the context of search problems where a team of searchers (traps) adopt the resetting 
strategy to make the search of the target more efficient [28]. For r = 0, it reduces to the 
standard target annihilation problem. For a nonzero r, the average survival probability 
in one dimension was computed exactly and it was found to decay as a power-law at 
late times [28] 

P 8 (t) « C rV 5 ^ (2) 

where C is a constant. For short times rt << 1, the decay is stretched-exponential [28] 
as in Eq. (JTJ). 

In this paper, we study analytically the average survival probability of a target 
particle in one dimension in another variant belonging to class (b). In this model, the 
target remains immobile at the origin as in the original target annihilation model. The 
traps also remain noninteracting and are distributed initially in space with a uniform 
density p. However, each trap, instead of performing a continuous-time Brownian motion 
with a diffusion constant D, now performs a continuous-time random walk (CTRW) 
with arbitrary jump length distribution and arbitrary waiting time distribution between 
successive jumps. Our result for the asymptotic survival probability includes the cases 
when the traps undergo subdiffusive as well as superdiffusive motion, thus generalizing 
the previously known result for subdiffusive traps [23] . 

In the standard version of CTRW, a trap moves by making successive random 
jumps at random times [29j [301 IS]- More precisely, starting from its initial position, a 
trap waits a certain amount of random time r (drawn from a normalized distribution 
ip{r)) and then jumps by a certain distance rj (also a random variable drawn from a 
continuous and symmetric distribution f(n)) (see Fig. (j2])). At its new position it again 
waits for a random time drawn from ip{r) (each time independently) and then jumps to 
another new position by a random distance drawn from f{ri) and the process continues. 
Thus, CTRW has two 'input' distributions, one for time and one for space, namely ipij) 
and f{rj). We also assume that these two distributions don't vary from trap to trap, all 
trap motions are characterized by the same ip(r) and f(rj). In this model, the immobile 
target particle at the origin gets annihilated whenever any trap trajectory goes past the 
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Figure 1 . Schematic example of a realization where the target at the origin is immobile 
and there are 6 noninteracting traps around it undergoing stochastic motion. The traps 
are initially placed uniformly in the box (— L/2, L/2). Eventually we are interested in 
the limit when the number of traps N — > oo, the box size L — > oo but with the density 
of traps p = N/L fixed. 



origin. Note that the particle at the origin gets annihilated even when a trap jumps 
over the origin. Thus the survival probability P s {t) of the target particle in this case is 
just the probability that not any of the traps jumps over the origin up to time t. Our 
goal is to understand how the two arbitrary input distributions ip(r) and f(rj) affect the 
canonical stretched exponential decay of P s (t) in Eq. (TjQ) that is valid when each trap 
performs continuous-time Brownian motion. 

The advantage of CTRW is that by appropriately choosing the two input 
distributions ipij) and f(rj), one can study various types of motions of a trap as special 
cases. For example, if one chooses ip(r) = S(r — At), one generates a discrete time jump 
process with jump length distribution /(r/) where the jumps happen instantaneously at 
time steps separated by At. If one further chooses a Gaussian jump length distribution 
f{rf) = exp [— r] 2 /2(r} 2 )} / ^j2-n(rf) and takes the limit (rf) — > 0, At — > but keeping the 
ratio (?7 2 )/(At) = 2D constant, one recovers the continuous-time Brownian motion. 
Similarly, by choosing a power law distribution f(rj) ~ M^" -1 for large \r]\ with 
< v < 2, one can generate Levy flights for a trap [30l ED [32] . Besides, the CTRW 
also provides a natural framework to generate anomalous diffusion process, namely 
subdiffusive and superdiffusive behavior that occurs in many natural processes [30j [31] . 
It would thus be interesting to know how the average survival probability P s {t) of a 
target particle decays with time when each of the traps, though still noninteracting, 
performs CTRW with arbitrary waiting time and jump length distributions, ip(r) and 
f(r)) respectively. 

Our main results are twofold: 

(i) We will first derive a very general result for the average survival probability P s {t) 
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for any model in one dimension belonging to class (b) above, i.e., when the target is 
immobile and each trap performs independently a stochastic motion Xi(t) starting 
from its initial position Xi(0). The stochastic evolution law of Xi(t) can be rather 
general. The only assumptions are (i) the process Xi(t) is invariant under a constant 
shift, i.e., Xi(t) and Xi(t) — c (where c is a constant) has the same evolution equation 
and (ii) the process Xi(t) —2^(0) is symmetric around the origin. Since the evolution 
law of the the z'-th trap Xi(t) — 2^(0) does not depend explicitly on the index i, we 
will generally denote it by Xt ra , p {t) = Xi(t) — Xi(0). For such a stochastic trap motion 
x tTap (t), we obtain a rather general result for the average survival probability 

P s (t) = exp[-2pE[M(t)}} (3) 

where E[M(t)] is the expected value of the maximum of the process x t v ap (r) over 
the time interval r £ [0, t], i.e, 

M(t) = max .[{xtrapCr)}] . (4) 

0<r<t 

We emphasize that this relation is very general and holds whatever be the stochastic 
process that describes the motion of a trap, as long as the process Xtm P (t) satisfies 
the two properties (i) and (ii) above. For example, if x tra _ p (t) is a continuous- 
time Brownian motion, it is well known that the probability distribution of 
the maximum M{t) of the process is a half-Gaussian [33] El]) i- e - 5 P(M,t) = 
exp[~M 2 /ADt}/VnDt for M > 0. Hence the expected value E[M{t)} = 2^/Dt/n. 
Plugging this result in the general relation in Eq. 03]), one recovers the classical 
stretched-exponential decay in Eq. (JT]). 

(ii) The relation in Eq. ([3]) holds, in particular, when x tra , p (t) represents a CTRW with 
arbitrary ip{r) an d f(v)- For such a general CTRW, we will first derive the exact 
asymptotic expression of E[M(t)} for large t. Plugging this result in the general 
relation ([3]) then provides us with the exact asymptotic behavior of P s (t) 

P 8 (t) « a exp[-M e ] (5) 

where the stretching exponent 6, the constant b and the amplitude a are determined 
exactly for arbitrary ip(r) and f{rf) (symmetric and continuous) (the precise results 
are detailed in Eq. ( )35|) ). We also verify numerically our analytical predictions. 

The rest of the paper is organized as follows. In section II, we derive the main 
results of this paper: section II. A provides the derivation of the general relation ([3]) 
while in section II. B, we derive E[M(t)] for CTRW by extending known results for 
E(M n ) in a discrete time setting [33J. In section III, these results are put together to 
obtain the asymptotic large t behavior of the survival probability P s (t) for a CTRW 
trap with arbitrary ip{r) and f{r]). Our general result will recover, as special cases, the 
known expression for traps performing diffusion and subdiffusion. In this section, we 
also compare our analytical results to numerical simulations and obtain good agreement. 
Finally we conclude in section IV with a summary and some open questions. 



Survival probability of an immobile target surrounded by mobile traps 
2. Derivation of the main results 



2.1. Survival probability for general stochastic motion of a trap 

We consider an immobile target at the origin in one dimension, surrounded by traps, 
initially distributed uniformly over space with density p and each performing an 
independent stochastic evolution Xi(t) starting from the initial position Xj(0). We will 
denote the process Xi(t) — X{(0) by a;t rap (t) that starts at the origin and satisfies two 
properties (i) and (ii) mentioned in the introduction. In order to derive an expression 
for the survival probability of the target at the origin, we follow the strategy used by 
Bray and Blythe [18] for the special case when x trap (t) is a Brownian motion around 
the origin. We consider first a finite number N of mobile traps starting at fixed initial 
positions {x(0)} = {xi(0), . . . , xtv(O)} on the interval (— L/2,L/2), cf. Fig. (JTJ and 
calculate the survival probability of the target up to time t for this fixed set of initial 
positions. Next we average this survival probability over the starting positions where 
each Xj(0) is drawn independently from a uniform measure on (— L/2, L/2) and take the 
limit N — > oo, L — > oo while keeping the density of traps N/L = p fixed. 

Let p(t\xi(0)) denote the probability that a single trap, starting at the initial 
position Xi(0), does not cross the origin up to time t. By independence of the N particles, 
the survival probability of the target factorizes into a product of the probabilities 
p(t\xi(0)) 

N 

Now averaging over all starting positions with uniform distribution on the interval 
(—L/2, L/2) and noting that particles were only distinguishable by their starting 
positions before averaging, 

N 

X\p{t\xM) 



s,N,L 



(t) 



=1 
1 

l -~L 



L/2 



dxp(t\x) 



L/2 



N 



L/2 



L/2 



dx (1 — p{t\x)) 



N 



(7) 



Taking the limit L — > oo, iV — > oo while keeping N/L 
p(t\x) = p{t\ — x) we get 



P s (t) = exp 

Note that, by definition. 
p(t\x) = Prob. 



2p / dx (l-p(t\x)) 



p fixed and using the symmetry 



(8) 



xAt) > for allO < r < t 



Xi{0) 



x 



(9) 



where Xj(r) represents the stochastic trajectory of a trap, say the i-th trap, starting 
at the initial position 2,(0) = x. Consider now the stochastic process x tr ap(T) = 
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Xiij) — Xj(0) = Xiij) — x. Then the probability p(t\x) in Eq. (J5J) translates into the 
following probability for the process £trap( T ) 

p{t\x) = Prob. £tra P (i~) > —x for all < r < t x trap (0) = (10) 

Using the spatial inversion symmetry of the process x tra p(r) around the origin, Eq. (ITU]) 
can be written as 

p{t\x) = Prob. x trap (r) < x for all < r < t x tTap (0) = (11) 

which, incidentally, is precisely the cumulative probability that the maximum M{t) = 
max < r < i [{x trap (r)}] of the process x trap (r) in the interval r G [0, t] is less than x, i.e., 

p(t\x) = Prob. [M(t) < x] . (12) 

Thus the probability density of the maximum is d x p(t\x) with x > 0. The expected value 
of the maximum is then given by E[M(t)} = J °° xd x p(t\x) dx. Noting that p{t\x) — > 1 
as x — )■ oo, it is actually useful to rewrite this as E[M(t)] = — f£°xd x [l — p{t\x)]dx 
and then do integration by parts giving 



E[M(t)] = / [l-p(t\x)]dx (13) 

Jo 

Substituting this result in Eq. ([8]) immediately gives the general result in Eq. (J3]). 

Let us make a quick remark here. For a random walk on a d-dimensional lattice, 

starting at the origin, the expected number of distinct sites visited by a walker up 

to time t is simply [35], ^^[1 — p(t\r)}, where p(t\r) is the probability that starting 

at r, the origin is not visited by the walker up to time t. For such a lattice 

walk in one dimension, the expected number of distinct sites visited up to t is thus 

S#[l — p(*I x )] = 2 XL>o[l ~ P(t\ x )]- Hence, the quantity ^ X>0 [1 — p(£|#)] ^ as ^ ne 
interpretation of half the average number of distinct sites visited by a walker up to time 
t. The lattice version of the relation ( ITBl) . where the integral is replaced by a sum over 
sites on the positive side, is consistent with this fact since in one dimension the average 
maximum up to t is precisely half the average number of distinct sites visited up to t, 
see 114, 361 and references therein. 



2.2. Expected Maximum E[M(t)} for the CTRW 

In this subsection, we focus on the case when x tra p(T) is a CTRW and derive the 
asymptotic large t expression for E[M(t)] of this process. The main idea behind this 
computation is as follows. Consider first a discrete-time random walk on a continuous 
line where the position x n of the walker at step n, starting at xq = 0, evolves at discrete 
integer time steps via the Markov rule 

X n = X n -x + T) n (14) 

where r] n 's are independent and identically distributed (i.i.d) random variables, each 
drawn a a symmetric and continuous probability distribution f{f]). Let M n = 
maxo<fc<n[{^fc}] denote the maximum of this process till step n. It turns out that 
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Figure 2. A typical trajectory xtrap(i") of a CTRW up to a total time t. The process 
starts at time r = at the origin and the first interval T\ also starts at r = 0. The 
trap waits at a point in space for a random time r drawn from a distribution ip( T ) 
and then jumps by a distance y drawn from a symmetric and continuous jump length 
distribution f{y). The successive time intervals {ri,T2,T n) T\ as t} spanning the total 
time t are statistically independent. 



an asymptotic expression for the expected maximum E[M n ] for large n can be derived 
explicitly for arbitrary symmetric and continuous jump length distribution f(rj) [33J. 
Knowing this explicit result for the discrete-time case, one can then derive the 
corresponding expression of E[M(t)] for large t for a CTRW by using a standard method 
involving renewal theory. In fact, it turns out that this renewal method, originally due 
to Montroll and Scher [29], is rather general in the sense that if one knows explicitly 
the expression for the expectation of any observable in the discrete-time jump process, 
the corresponding expectation value for the CTRW can be derived in a straightforward 
manner, as illustrated below. This fact sometimes goes by the name 'subordination 
property' in the literature [31J. 



2.2.1. Subordination property: Let O n be the expected value of any observable 
associated with the discrete-time process in Eq. ( THl) . For a corresponding CTRW, 
these steps do not occur 'at every tick of the clock', but rather at random points in 
time separated by random variables {t„}, called 'waiting times'. The waiting times r n 's 
are i.i.d variables each drawn from a waiting time density (WTD) denoted by ip(r). 
For a fixed total time t, the number of 'kicks' n is therefore a random variable and let 
Q(n\t) denote the probability that there are exactly n 'kicks' in time t. The continuous- 
time expectation value 0(t) of the observable in CTRW is then related to its discrete 
counterpart O n by the relation 



0(t) = J20 n Q(n\t). 



(15) 



?1=0 
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To compute Q{n\t) one can use the renewal property, i.e., the fact that successive 
intervals are statistically independent. Consider the interval [0, £], broken up into 
(n + 1) intervals (see Fig. (J2J)). The first n intervals are denoted by Ti, T2, . . ., r n 
and the last interval ri ast = t — ^27=1 Ti which is f ree °f 'kicks'. While each of Tj's 
(for i = 1,2, ...,n) is drawn from ip(r), the distribution of the last interval ri ast is 
slightly different (see Fig. (J5J)). The probability that an interval r is 'kick' free is simply 
4>(r) = f°° ip(r) dr = 1 — L i/j(t') dr' . Thus, ri ast is drawn from the density function 
4>(t), while the preceding ones are each drawn independently from ip{ T )- Then the 
probability Q(n\t) is given by the convolution 



Q(n\t) = ... I dridr 2 ...dr n dri ast ip(n) t/i(t 2 ) • • • ip{r n ) </>(riast) 

Jo Jo 

x 6{ti + r 2 + ... + r n + r^t - t) . (16) 

Convolutions factorize under a Laplace transform. Denoting the Laplace transform of 
the WTD by ip(u) and noting that J °° dt e~ ut <p(t) = (1 — ip(u))/u, one gets 

0(u)= r dt e~ ut Oft) = 1 ~ ^ U > Y* O n U(u)] U . (17) 

Jo « fo L J 

Now since < ip < I for u > 0, the sum on the right hand side of Eq. (j!7p can 
be interpreted as the generating function h(s) = Y^=o^nS n of O n with s = ^(u). 
Then the Laplace transform 0(u) of the continuous-time quantity 0(t) is related to the 
generating function h(s) of the discrete-time process via Eq. (ffj 



2.2.2. Typical displacement of a trap: As an example of the subordination method, let 
us consider a discrete-time process ( Tl4l) where the jump length distribution f(rj) has a 
power- law tail for large \r)\ 



m 



1^11+!/ 



with v > 0. The variance of the jump length distribution a 2 = J rf f(rj) drj is finite if 
v > 2 and infinite if v < 2. For z/ < 2, this process is the so called Levy flights [52" } I5U I I3T] . 
Let P n (x) denote the distribution of the position of the trap after n steps, starting at 
the origin. It is then well known and is easy to show that P n (x), for large n, has the 
scaling form [30] 



*M = ;^U0 (19) 

where the scaling function, for large \z\, decays as C^z) ~ |^| _1_M for < /i < 2 and is 
Gaussian for jjl > 2. Thus, the typical displacement of the trap after step n scales, for 
large n, as 

a; typ( ri ) ~ nl f° r < v < 2 (superdiffusive) (20) 



r^j 



n 1/2 for v > 2 (diffusive) (21) 
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Let us also consider a WTD ip(r) that has a power-law tail for large r 



-l+a 



(22) 



with a > 0. The mean waiting time fi = L Tip(r) dr is finite if a > 1 and infinite if 
a < 1. Correspondingly, the leading asymptotic behavior of the Laplace transform V>( M ) 
of the WTD can be computed near u — > (see, e.g., appendix 2 of Ref. [37]) as 

{1 — n u for finite /i 
(23) 
I — A \T(—a)\u a for infinite /i 

Replacing O n by a; typ (n) in (fTTj) and using the above asymptotic expressions one can 
easily check that x typ (t), the typical displacement of the trap undergoing CTRW, behaves 
asymptotically in different ways depending on the following four cases 

(i) both a and /i finite 

(ii) a finite, /i infinite 

(iii) a infinite but \x finite 

(iv) both a and \x infinite 

which, for power- law jump length distribution j{r\) in Eq. ( TT5|) and power-law WTD in 
Eq. (|22l) correspond to 

(i) v > 2 and a > 1 

(ii) v > 2 and a < 1 

(iii) z/ < 2 and a > 1 

(iv) v < 2 and a < 1. 

This gives 

t 1 / 2 case 1 



x typ {t) ~ < 



t a/2 case 2 
t 1/!/ case 3 



(24) 



t a/v case 4 

where the cases refer to the four different situations mentioned above. Note, in 
particular, that by appropriately choosing the input exponents a and u, one can generate 
diffusive (case 1), subdiffusive (case 2), super diffusive (case 3), as well as a completely 
anomalous (case 4) where one can get both subdiffusive and super diffusive behavior by 
appropriately choosing a/ v. 

As an aside, another recent example where this subordination property has been 
used concerns the statistics of records. For a discrete-time jump process in Eq. (IT4"|) . 
the statistics of the number of records up to step n turn out to be completely universal, 
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i.e., independent of the jump length distribution f(rj) as long it is symmetric and 
continuous [38] . For instance, the average number of records up to step n, grows as 
yln/V for large n J3S]- Using the subordination property mentioned above and 



(Rn) 

by choosing O n 



R n , the average number of records up to time t for a CTRW with 



power-law WTD as in Eq. ( 1221) has been computed recently 



2.2.3. Expected maximum of the trap motion: Having obtained the late time scaling of 
the typical displacement of a trap undergoing CTRW, we now turn to computing the 
expected maximum E[M(t)} of the trap motion Xt T &p{t)- In order to use the general 
subordination technique mentioned above, we then need to know the behavior of E[M n ] 
for large n of the discrete-time process (|14p where M n denote the maximum up to step 
n. Fortunately this can be done as illustrated below. 

Computing explicitly the full distribution of M n for arbitrary f(jj) is, in general, 
a hard problem [33l [31]. However, there exists a general formula known as Pollaczek- 
Spitzer formula [301 EH] that reads [33] 



J2s n E[ 



-pM n ] 



n=0 



Vl^ 



<p{s,p) 



where 



(f)(s, p) = exp 



P_ 

n 



In (1 - sf(k) 

p 2 + k 2 



dk 



(25) 



(26) 



and f(k) = j_ f(rj) e lkv d-q is the Fourier transform of f{i]). From the expression 

the generating function of the expected maximum can, in principle, be obtained via the 

formula 

1 d(f)(s,p) 



h(s) =^s n E[M n 



n=0 



p=0 



VI - s dp 

Replacing O n by -E[M n ] in the general formula ( IT71) . we then have for CTRW 

l-i>{u) 



dte- ut E[M(t)} 



u 



h[ip(u] 



(27) 



(28) 



To extract the leading asymptotic behavior for large t, one needs to analyze the 
right hand side of Eq. ( 1281) in the limit u — > 0, or equivalently h(s) in Eq. ( 1271) in the 
limit s — > 1 . Extracting the leading singularity of the right hand side of (|27|) in the limit 
s — > 1 turns out to be rather nontrivial [33]. For jump length distributions with finite 
second moment a 2 = J rf f(rj) drj, the principal singular behavior of h(s) near s — 1 
can be extracted explicitly 



h(s) 



a 



1 



72(l-s) 3 / 2 tt(1-s 



dk 



In 



/(*) 



a*- 
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In case when a 2 is infinite, as for instance the case for Levy flights in Eq. ( fl8|) with 
v < 2, where f(k) w 1 — \ak\ v to leading order for k — > 0, the asymptotic behavior of 
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Figure 3. The four different asymptotic power-law growth of E[M(t)] for large t in 
the four cases in the (a — v) plane is summarized. 



h(s) near s = 1 turns out to be different from Eq. ([291 . For example, for 1 < v < 2, 
one can show [33] 
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where B(a, b) = T[a]T[b]/T[a + b] is the standard Beta function. 

For the case < v < 1, it turns out £[M n ] is infinite for any n > 0. Subsequently 
i~7[M(£)] is also infinite for any finite t. Consequently, for < v < 1, the survival 
probability P s (t), via the relation ([3]), is for any finite time t (the target is killed 
instantaneously). This rather pathological case will not be discussed further. 

We now need to substitute these asymptotic s — y 1 behaviors (for the two cases 
in Eqs. ((29]) and (|30|)) in Eq. ([28]) and analyse the large t behavior of E[M(t)]. To 
carry out the asymptotic analysis, one again needs to distinguish four different cases 
mentioned in the previous subsection, namely (1) both a and \x finite (2) a finite but 
fj, infinite (3) a infinite but \x finite (4) both a and /x infinite, where a 2 and ji refer 
respectively to the second moment of the jump length distribution f(n) and the first 
moment of the WTD ip(r). When both f(rj) and ip{r) have power-law tails as in Eqs. 
(II 8p and ([22]) respectively, the above 4 cases correspond to 4 regions in the (a — v) plane 
(1) v > 2 and a > 1 (2) v > 2, a < 1 (3) v < 2, a > 1 and (4) v < 2, a < 1. As 
mentioned before, we will not discuss the pathological case < v < 1 and will restrict 
ourselves in the (a — v) plane only for v > 1 and a > (see Fig. ([3])). 

The next step is to plug the appropriate asymptotic expression for the Laplace 
transform of the WTD from Eq. ([23]) into the corresponding expression for the 
generating function of E(M n ) in Eqs. (|29l and ([30]) and subsequently invert the Laplace 
transform in the asymptotic large t limit. One needs to use the Tauberian theorem, cf. 
[4*2] [43] that states that a Laplace space behavior of u~ p for m< 1 corresponds to a t p_1 
behavior for i > 1 under inversion of the Laplace transform. It turns out that three 



(30) 
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constants appear in the final expressions for E[M(t)] which we denote below as 
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with B(a, b) being the standard Beta function. In the definition of Eq. ( |3T1) we assume 
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v > 2 (a 2 finite) and in Eq. (I3"2"j) . one assumes 1 < v < 2. We then obtain for large t 
the following asymptotic expressions for E[M(t)} in the four cases 

^k t 1 ^ + # case 1 
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where the cases refer to the four different situations mentioned above. Note that the 
result in case 1 also holds for arbitrary distributions f(j]) and ip{r) such that a 2 and 
jjl are finite. These 4 different asymptotic power- laws for the growth of E[M(t)} in the 
(a — v) plane are summarized in Fig. ([3]). Note that the leading asymptotic growth 
of E[M(t)} for large t is similar to that of the typical displacement in Eq. (|2"4"j) in all 
the 4 regions, except the prefactors of the power-law growth of E[M(t)} and x typ (t) are 
different. As in the typical case, E[M(t)} shows the usual diffusive scaling to leading 
order in case 1, subdiffusive in case 2, super diffusive in case 3, while in case 4 it scales 
like t a l v . Interestingly, this totally anomalous case 4 of Levy flights interspersed with 
long waiting times can exhibit diffusive scaling, for any choice of a and 1 < v < 2 but 
with ajv = 1/2. We have verified the analytical predictions in Eq. ( 1341) via numerical 
simulations. The results are shown in fig H] where different types of anomalous scaling 
behaviors are compared. Numerical results are in excellent agreement with analytical 
predictions. 



3. Asymptotic Survival Probability: Analytical Predictions and Numerical 
Simulations 

Plugging the expressions for E[M(t)] for the four cases considered above in the general 
relation in Eq. ([3]), we obtain the following large t behavior for the average survival 
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Figure 4. Expected maximum of a CTRW as function of t. On the log-log scale, the 
exponent of t can be read off from the slope. Fig. (4a) corresponds to the case of finite 
/i while Fig. (4b) corresponds to infinite fj,. The bottom curve of Fig. (4a) corresponds 
to case 1 where both /j, and a are finite and the exponent is 1/2 (see Eq. (|34]l ). The 
top curve of Fig. (4a) corresponds to case 3 where /i is finite but a is infinite- the 
exponent is 1/v. Similarly, in Fig. (4b) the bottom curve corresponds to case 2 (where 
(i is infinite but a is finite and the exponent is a/2 as in Eq. (|34[l ) while the top curve 
corresponds to case 4 (where both fi and a are infinite and the exponent is a/v). Note, 
in particular, that in the bottom curve of Fig. (4b), by choosing a/v — 1/2 one obtains 
the diffusive behavior even though the process is totally anomalous. Simulations used 
a Levy-stable or standard normal density for the jump lengths and either exponential 
or Pareto distributed waiting times. Parameters are given in the figure. 
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(35) 



exp [-2p {Jt a ' u + I)] case 4 

where the constants K, I and J are given respectively in Eqs. (I3"T|) . (|33|) and (1331 . Thus, 
in all four cases, one finds generically the stretched exponential behavior for large t 

P s (t)^aexp[-bt e ] (36) 

except that the amplitude a, the constant b and the stretching exponent 9 differs in the 
4 cases. In the 4 different cases, the exponent 6 is respectively (1) 9 = 1/2 (2) 9 = a/2 
with a < 1 (3) 9 — 1/V with 1 < v < 2 and (4) 9 = a/v with a < 1 and 1 < u < 2. 
We also emphasize that our method yields an exact result even for the amplitude a in 
front of the stretched exponential in Eq. ( |36j) . Usually this is much harder to obtain 
as it corresponds to computing the subleading correction term in E[M(t)}. In the four 
respective cases, we get the amplitudes 
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(i) a = exp[— 2 pK] 

(ii) a = exp[— 2 pK] 

(iii) a = exp[— 2 pi] 

(iv) a = exp[— 2 pi] 

where the constants K and / are given respectively in Eqs. (J3"T|) and fl32|) . 
3.1. Recovering known results for diffusive and subdiffusive traps 



From our general result in Eq. ( 1351) . it is easy to recover the known results for the 
Brownian (diffusive) and subdiffusive traps first derived in [Hj and [23] respectively. 
Consider first the Brownian case. To recover this result corresponding to Brownian 
motion of a trap, we need to focus on case 1 of our result in Eq. ( 135]) and take the 
Brownian scaling limit where both mean square jump length a 2 and mean waiting time 
p tend to zero, while their ratio D = o 2 j2p stays constant. In order to take this limit 
cleanly, it is convenient to first make a change of variable k — y ak/\/2 in the expression 
for the constant K in ( 13T1) to get 

Note that K is independent of p. Since f(k) — y as k — y oo, we can safely take the 
limit a — y in Eq. (|37|) and obtain K = 0. From case 1 of Eq. (|35|) . we then recover 
the Brownian result: P s (t) —y exp(—2p^ADt/7r). 

To recover the result for the continuous-time subdiffusive case in [23], we need to 
take the scaling limit of our result in Eq. ( !3~5l) corresponding to case 2. As in the 
Brownian case, we need to take the a 2 — y limit. However, unlike in the Brownian 
case, here the mean waiting time p is already infinite. So, the only way a sensible limit 
of line 2 in Eq. (13 5 j) can be reached if 2y4|T(— a)| goes to zero while a 2 — y 0, but their 
ratio D a = a 2 /(2A\T(— a)\) remains fixed. This ratio D a can be called the generalized 
diffusion constant since it can be shown that it appears also in the amplitude of the 
mean-square displacement in case 2 (with a < 1) 

(4,(0) -F^f. (38) 

This relation can be thought of as the generalization of the standard diffusion relation 
( x trap(t)) — 2Dt with D being the standard diffusion constant. Indeed Eq. (|38l) reduces 
to the standard diffusion relation when a — y 1. Taking this scaling limit in line 2 of Eq. 
(135]) , we thus recover the result of Yuste and Lindenberg [23] 



P s (t) ->• exp 



2py/lX a/2 
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(39) 



r(l + a/2) 

Our general result in Eq. (13"5"]) . apart from recovering known results in the 
appropriate limits as shown above, also provides new results for the super diffusive case 
(which can be achieved in cases 3 and 4). In addition, our result also provides not just 
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Figure 5. Left: For standard normal distributed jumps length with a 2 = 1. This 
corresponds to case 1. The dotted line shows the analytical behavior Right: For Levy 
jump length distribution with Fourier transform f(k) — exp[— (a|fc|) y ] with the choice 
ix = 1.5 and a = 1. While for the Gaussian case, the survival probability was averaged 
over 10 5 realizations, for the Levy-distribution, between 10 6 (AT = 10, L = 1000) and 
10 7 (all others) realizations were used, since due to the possibly large jumps, only very 
little realizations survive up to long times. For the same reason, the density of traps 
is a factor of 10 lower than for the Gaussian case. 



the leading asymptotic behaviors in all cases, but also the subleading order corrections 
in the form of non-universal amplitudes. 



3.2. Numerical simulations 

To confirm our analytical results, we also performed numerical simulations. Due to the 
computational effort of simulating large systems, only a finite number of iV particles 
uniformly distributed over the interval (0, L) with N/L = p fixed were considered. 
The reason for only considering traps on one side of the origin in just computational 
convenience. For better accuracy of comparison, fig [5] compares — ln(P s (t))/t e to 
numerical simulations with the exponent 9 taking the value corresponding to the case 
considered, see Eq. (|35|) . In all cases, — \ia(P s (t))/t e is of the form pa x + pa 2 /t e with 
constants a± and a 2 depending on the cases. The missing prefactor 2 infront of the trap 
density p is due to the asymmetry of having only traps on one side of the origin. 

Simulations were then performed for increasing values of iV and L to show that as 
the effects of finite system size become less important, the simulated curves converge 
onto the analytical predictions, see Fig. [5j The simulations were performed with the 
WTD ip(r) — S(r — 1), thus leaving us with only cases 1 and 3 to distinguish. The 
influence of infinite mean waiting time on the expected value of the maximum of a 
CTRW has already been illustrated numerically in Fig. HI so it suffices to study the 
survival probability of a target surrounded by traps performing a random walk or Levy 
flight with finite p. From Fig. [51 one sees that with increasing size L, the numerical 
curves approach the analytically predicted behaviors shown by dotted lines. 
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4. Conclusion 

To summarize, we have studied analytically, in one dimension, the survival probability 
P s (t) of an immobile target at the origin surrounded by a sea of independent mobile 
traps initially distributed uniformly with density p on the line. We have shown that for 
a rather general stochastic motion Xt rap (t) of each trap, there is a very simple general 
relation, P s {t) = exp [— 2pE[M(t)]], that relates the survival probability to the expected 
maximum of the process Xtrap{t)- This general result recovers all earlier known special 
cases of x tTa _ p (t). In addition, we have presented new result when x trap (t) is a continuous- 
time random walk (CTRW) with arbitrary jump length distribution f(r]) and with 
arbitrary waiting time distribution ip(r). By choosing these two input distributions of 
CTRW one can generate a variety of stochastic motions of the trap, both subdiffusive 
and super diffusive. For such a general CTRW, we were able to evaluate exactly the two 
leading terms in the asymptotic behavior of E[M(t)} for large t and plugging this in our 
general relation, we were able to compute the two leading terms for the large t behavior 
of P s (t). Generically, we found that for large t 

P a (t) ~a exp[-bt e ] (40) 

where the exponent 8, as well as the two constants a and b are computed exactly in this 
paper when x trap (t) is a CTRW with arbitrary f(r]) and ip(r). Previous results in the 
literature were known only for diffusive and subdiffusive x tvap (t) (case 1 and case 2 in 
Fig. (J3])). Our result is more general and in particular, the exact asymptotic expressions 
for the two remaining cases (the super diffusive (case 3) and totally anomalous (case 4)) 
are, to our knowledge, new. In addition, our method allows us to compute even the 
nonuniversal amplitude a (and not just the stretching exponent 9) exactly in all cases. 

There are several future directions in which it may be possible to extend our results. 
First, it would be interesting to generalize our result to higher dimensions. The main 
difficulty in higher dimensions is that it is difficult to obtain an explicit expression for 
the probability that a discrete-time random walker does not hit a static ball (target) of 
finite radius up to time £ [9] , though it may be possible to obtain some exact results 
in the special case of d = 3 [H QUI QTJ . It would be interesting to derive the result for 
P s (t) in this more realistic d = 3 case. 

Secondly, it would be interesting to derive the corresponding result for P s (t), in 
one or higher dimensions, in the case when the target itself is mobile, e.g., when the 
target also performs a general CTRW as the trap. As mentioned in the introduction, 
in special cases of diffusive and subdiffusive motion of the target it is known that the 
leading asymptotic behavior (but not the subleading) of P s (t) does not change from the 
situation when the target is static. Is this fact still true when the target performs a 
general CTRW? It would be interesting to prove or disprove this fact. 

Finally, it would be interesting to find suitable applications of our results. In fact, 
the totally anomalous case 4 (see Fig. ([3])) for which we present new exact results, may 
at first glance seem to be only of academic interest, but has recently been realized in 
a more physical context, namely, in a model for spread of epidemics by human travel 
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J. The results presented above thus have a natural interpretation as the probability 
of staying on the 'safe' side of the propagating front of a disease up to time t and it 
would be interesting to explore in detail this interesting possible application in future. 
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